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INTRODUCTION 


This  work  is  aimed  at  the  understanding,  description,  and 
prediction  of  magnetohydro  dynamic  phenomena  exhibited  under  conditions 
of  extremely  high  interaction  and  large  magnetic  Reynolds  number.  The 
plasmas  which  consistute  the  fluid  medium  in  such  flows  may  exhibit 
nonideal  thermodynamic- and  kinetic  behavior.  The  theoretical  work  at 
-,STD~-Reaeaxch- under  QNIL-aupport  has  two  principal  objectives:  ft )  the 
elucidation  of  basic  phenomena  in  strong  interaction  high  magnetic 
Reynolds  number  flows  independently  of  specific  experiments  or  machines; 
and -fift  the  perfection  of  predictive  theories  to  accurately  describe  and 
model  specific  experiments  aimed  at  magnetohydrodynamic  power 
production. 


In  what  follows  we  present  the  general  mathematical  description 
of  magneto  gas  dynamic  flows  in  the  high  magnetic  Reynolds  number  regime. 
We  present  several  illustrative  calculations  of  quasi-one -dimensional 
transient  effects  in  strong  interaction  flows.  -Wethen  pro  sent  two- 
dimensional  high  Reynolds  number  electricity  results^  including  the 
realistic  effects  of  nonuniform  velocity  and  electrical  conductivity 
resulting  from  hypersonic  boundary  layers  and  from  shock-induced 
nonunifo  r  mitie  s .  . 


2.  THE  MATHEMATICAL  DESCRIPTION  OF  MAGNET OGASDYNAMIC 

FLOWS 

2.  i  Fluid  Conservation  Laws 

We  may  describe  the  fluid  in  terms  of  its  mass  density,  p, 
velocity  O’,  and  internal  energy  e.  We  describe  the  electromagnetic 
effects  in  terms  of  the  electric  field  E  and  magnetic  field  B*.  These 
variables  are  considered  to  be  general  functions  of  space  x  and  time  t. 

The  conservation  laws  for  mass,  momentum,  and  energy  are 


*  (pU)  =  0  (i) 

(pU)  +  (pUU)  =  7-T  +  JxB  (2) 

^  [p(e  +  U2/2)]  +  V-[(e  +U2/2)  Pu]  *  V*  (¥•  u)  -V-q+7  •  E  (3) 


In  the  conservation  laws,  7T  is  the  total  pressure  tensor  and 
q  is  the  heat  flux  vector.  This  system  of  conservation  laws  is  completed 
in  the  limit  of  infinitely  fast  kinetics  by  the  kinetic  and  caloric  equations 
of  state 


P  =  p(p  ,  «  ) 

(4) 

e  =  e  (P.T) 

(5) 

where  p  is  the  isotropic  part  of  the  stress  tensor  IT  and  T  is  the 
temperature.  For  a  general  fluid  the  state  equations,  Eqs.  (4),  (5)  cannot  be 
explicitly  given  but  are  embedded  in  the  general  statistical  mechanical 
description  of  the  equilibrium  thermochemistry  of  the  system. 
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In  the  case  of  a  perfect  gas  with  particular  gas  constant  R 
and  specific  heat  ratio  y  explicit  formulae  may  be  given: 


p  =  pR  T 

(6) 

e  =  ( y-1)"1  R  T 

(7) 

p  -  p(yOe 

(8) 

2.2  The  Electromagnetic  Contributions 

The  electrical  equations  (consisting  of  the  Maxwell  equations 
and  the  generalized  Ohm's  law)  govern  the  electric  and  magnetic  fields 
£,B  and  the  conduction  current  density  7.  In  the  hydromagnetic  limit 
these  are 


rj  y  TT  _  9B 

VXE  -  "  9t 

(9) 

VXB  =  h-qJ 

(10) 

V  •  B  =0 

(ID 

7  =  <t(e+uxb)+7k 

(12) 

where 

is  the  thermal  diffusion  flux  vector, 

7K=  aK 

and  K 

is  given  by  1 1  ] 

■ 

K  s  - 

9(1)\7T  +e*2)VT  XB  +9(3)(VT  XB)XB 
e  e  e 

N 

•  E  r  )XB 

or=l  L  J 

(13) 
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f.  ! 


The  ••  •  and  |3^  •  •  •  are  transport  coefficients  defined  in 
[l],  Tg  is  the  electron  temperature  and  the  subscript  a  denotes  a 
plasma  component. 

The  contributions  to  the  momentum  and  energy  of  the  system 
by  the  electromagnetic  field  are  contained  within  the  Lorentz  force, 

J  X  B  and  the  Lorentz  power  J  •  E  .  The  Lorentz  force  may  be 
represented  in  terms  of  the  Maxwell  stress  tensor  T  as 

7  XB  =  7-  T  (U 

where  the  Maxwell  stress  tensor  is  defined  as 

7  =  £q‘  b  -  bz/2  T)  df 

Correspondingly  the  Lorentz  power  may  be  represented  in  terms  of  the 
Poynting  flux  S  and  the  electromagnetic  energy  density  e^ 

J  .E  =  -V  •  r  -  (1< 

The  Poynting  flux  5  is  defined  as 

7  =  Cq"  1  E  X  B  (11 

while  the  electromagnetic  energy  density  is 

=  ~o‘  b2/2  1“ 

Let  us  expand  the  Poynting  flux  in  terms  of  J  ,  B  through  the  use  of 
the  Ohm’s  law 


=  <r( 


E  +  U  X  B  +  K 


We  have 
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~S  =  ^‘(exb)  =  (p.Q<rJ  JxB  -^Q1  (uxb)xB  -^KXB 

The  term  JXB  is  simply  V  •  ~  .  The  term  (U  X  £*)  X  B  is  readily 
shown  to  be 

(UXB)XB  =  U  •  (BB)-D(B'B) 
which  can  be  rearranged  to 


(UXB)XB  =  U  •  (B  B  -  B2/2  I  )  -  U  (B2/2) 


The  Poynting  flux  is  therefore  represented  as 


s  =  e  U-U'f+piJ-r-nJ^XB 

m  j\ 


(19) 


/V  -  1 

where  r)  =  (jx^ <r)  is  the  magnetic  diffusivity.  We  note  that  the  Poynting 
flux  can  be  decomposed  into  four  constituent  parts:  (a)  a  purely  convected 
flux  of  electromagnetic  energy  carried  by  the  motion  of  the  medium 
(emU );  (b)  a  power  flow  represented  by  work  done  per  unit  time  by  the 
Maxwell  stresses  acting  on  the  moving  medium  -U  mT  ;  (c)  a  diffusive 

flux  of  electromagnetic  energy  driven  by  gradients  of  the  Maxwell  stress 
tensor  and  (d)  a  power  flow  represented  by  work  done  per  unit 

time  by  the  thermal  diffusion  flux  ("OJ^  *  B). 

The  electric  and  magnetic  fields  and  currents  may  be  expressed 
in  terms  of  vecter  potentials  A  and  scalar  potential  $  as 


B  =  V  x  A 


E  =  -V  <5 


9  A 
9t 


(20a) 

(20b) 
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2.3 


viscous , 

The  mass,  momentum,  and  energy  equations  for  the 
hydromagnetic  system  may  now  be  expressed  as 

general, 

Iftv-sr  =  o 

(21) 

Qjl  +  v  -  r  =  v  •  V 
at 

(22) 

•  IT  =  V  •  (u  •  F)  -  V  •  (r,V.  T)  -  V'q  +  V-ln^XB 

(23) 

In  the  above  M  -  pU  is  the  momentum  density  and  T  is  the  total  fluid 
and  electromagnetic  momentum  flux 

F  =  pUlf  +  pf  -  F  (24) 


The  total  energy  density  is 

e  =  p(e  +  U2/2)  +  e  m  (25) 

and  s  is  the  total  enthalpy  flux  vector 


The  fluid  stress  tensor  it  has  been  decomposed  into  a  pressure  p  and 
viscous  stress  tensor  r  where 


P 


j  Trace 


(ff  ) 


(27) 


and 


it 


+  T 


(28) 


We  may  define  the  electromagnetic  pressure  pm  as  the  mean 
normal  compressive  Maxwell  stress: 


6 


(29) 


p  =  "  T  Trace  (  T  ) 
m  5 


Expressed  in  terms  of  the  magnetic  field  intensity  B  ,  the  pressure  is 

Pm  =  3(B2/2^o)  (30) 


or  expressed  in  terms  of  the  electromagnetic  energy  density  e 


m 


Pm  ~  3 


m 


(31) 


The  Maxwell  stress  tensor  may  be  decomposed  into  a  mag¬ 
netic  pressure  p  and  a  Maxwell  stress  deviator  from  isotropy  T *  as 
r  m 


T  ~  T*  ~  Pm  1 


(32) 


Equation  (32)  may  be  thought  of  as  the  defining  equation  for  the  Maxwell 
stress  deviators  Y*.  The  mass,  momentum,  and  energy  equations  may 
now  be  rewritten  in  alternative  form  as 

||  +  V  •  M  =0  (33) 

V-7  (34) 

^  +  7  •  H  =  7  •  (U *  1?*)  +  V  '(U-*rj  -  V  •q’  -  V  -  (nV  •  F)  +  7  •  (r)7KXB  )  (35) 

In  the  above,  G  is  the  total  momentum  flux  with  only  the  magnetic  pressure 

included  as  the  magnetic  contribution 

G  =  pUU  +  (p  +  pm)7  (36) 

while  H  is  the  total  enthalpy  flux  consisting  of  both  fluid  and  electro¬ 
magnetic  pressure  and  energy  contributions. 

H  *  [e  ♦  (p  ♦  pm)]  U  (37) 
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In  terms  of  the  total  energy  density  e,  the  momentum  density 
M ,  and  the  electromagnetic  energy  density  the  state  equations  (6)  - 

(8)  become 

p  =  (Y-l)[e  -  em  -tf2/2p]  (38) 

T  =  p/pR  =  £e  -  em  -  A/2/ 2pJ  (39) 

Let  us  now  consider  the  transformation  of  the  electrical 
equations  (9) - (i 2)  into  more  useful  forms.  Combining  Eqs.  (9)-(12)  we 
obtain  the  governing  equation  for  the  magnetic  induction  B: 

^  -  V  x  (U  XB)  =  -7  x  (rjV  X  if)  +  V  x  K  (40) 

We  note  that  given  the  magnetic  induction  B(x,t)  governed  by  Eq.  (40) 
one  immediately  has  specified  the  Maxwell  stress  tensor  T  and  the 
electromagnetic  energy  density  e^.  Further,  the  current  density  J  is 
determined  from  B  as 

7  =  £0~iv  X  B 

and  the  electric  field  E  as 

E  =  -U  X  B  +  r,V  X  B  -  K 

2. 4  Viscous  and  Heat  Conduction  Effects 

Let  us  now  make  some  observations  about  the  viscous  stress 
tensor  t  and  the  heat  flux  vector  q  .  The  Navier-Stokes  moments  of 
the  Boltzmann  equation  yield  kinetic  theory  forms  of  these  quantities: 


tL  = 

(41) 

■q,  =  -WT 

(42) 
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In  the  above,  p.  and ^  X.  are  the  coefficients  of  viscosity  and  thermal 

conductivity  and  (VU)Q  is  the  symmetrized,  traceless  velocity  gradient 
tensor: 


/3U.  3U:, 


3x, 


ij 


We  denote  the  stress  tensor  and  heat  flux  vector  with  a  subscript  L  to 
denote  that  these  are  laminar  quantities.  If,  on  the  other  hand,  we 
interpret  the  fluid  variables  p,U,e,T,  as  turbulent  mean  quantities, 
then  t  and  q  contain  turbulent  contributions  due  to  the  turbulent  velocity 
and  enthalpy  correlations.  Hence,  the  complete  stress  and  heat  flux 
fields  for  a  turbulent  hydromagnetic  medium  are 


7  =  <Pu'  u')  +  7l 

«T  =  <  pU  h  )  +  ~qR  +  Tl 

— r 

where  U  ,h  are  the  turbulently  fluctuating  velocity  and  enthalpy  and  <> 
denotes  an  ensemble  average.  A  detailed  higher  order  closure  theory  for 
the  turbulent  contributions  (pU  U  ),  (  pU  h  )  is  given  by  Demetriades, 
Argyropoulos,  and  Lackner  [2].  The  radiative  heat  flux  is  q^.  Since 
the  optical  depths  in  dense,  explosion  generated  plasma  are  so  smali, 
the  radiative  heat  flux  is  only  important  in  layers  near  the  plasma  surface 
of  the  order  of  the  radiation  free  path. 

2.  5  Nondime  ns  ional  Forms 

Let  us  consider  the  nondimens ionalization  of  the  fluid  equations 
(33)  -  (35)  and  the  electrical  equation  (4  ).  For  this  purpose  let  us 
specify  characteristic  values  of  the  variables  as  pQ, Uq.Cq,  pQ,  TQ,  •  •  •  as 
well  as  magnetic  field  BQ.  We  define  a  characteristic  length  L  and 
characteristic  time  tQ  =  L/Ufl.  We  indicate  nondimensional  variables 
with  f). 
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The  fluid  and  electrical  conservation  laws  then  become 


^  +  V  •  M 


0 


__  at 

+  v  •  f  =  R"1  v  • 
at  e 

2M7U)0J 

-1 

%  +  V  •  S  =  R’ 1  V  •  F  2p  U  •  (V  U)0 

at  e  L 

]+[(Y-l)M2RePR]  V 

-  SR^v-(nV-f)  - 

2SV  .(rjJgXB  > 

^  ^  <v 

•  (\  V  T)  - 


(43) 

(44) 


(45) 


—■  •  V  X  (U  X  B)  :  -R'1  IvxfiVXB)  -  7  x  ( n  J„) 

at  K 


B  =  7  x  A 


at 


(46) 


In  the  form  in  which  the  Maxwell  stress  tensor  is  decomposed  the  left  hand 
sides  of  the  momentum  and  energy  equations  take  the  forms 


7K  U  ^  ^  «w«w  CX  m 

+  v  .  g  -sv  •  r*  =  “  +  v  •  r 
at  at 


~  +  V  •  H  -  SV  •  ( U  •  T  )  =  +  V  •  ff 

at  at 


The  nondimensional  variables  are  defined  as 


t  =  U0t//L  ^  =  L7 

U  =  u/u0 

P  =  p/pQ 

5  =  «/p0U0 
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r 

.  t 

C 

[  T  =  T/T„ 

]  e  =  p  je  [v(‘^)m2]  +  U2  j-  +  S  em 

j  *  €  =  e/tW-lf1^] 

:  <[  «m  =  (B2/2?0)/(b2/2{20) 

!  E  =  E/U0B0 

•  B  =  B/B 0 

1  A  =  A/(Bq/L) 

*  =  «/(U0B0L) 

J  =7/(^1B0/L) 

=  JK^°0U0B0* 

G  =  pU  U  +  j^(\M2)  p  +  SpJ  T 
f  a  p  U  U  +  hM2)*1  pT  -ST 

*  fm/(Bo/2^) 

if  .  [;  +  (vm2)  p  +  s;J  5 

a  s  [(e  +  (\M2)  p)T  -  s  T  j  -  U 
T  =  r/(B2/2p0) 

v  =  »*  =  k/V0  n  =  n/r)Q 

where 

X0  =  X(p0’T0)’  »»o  *|i(pofV»  no  =n(po,To)* 

The  state  equations  (38),  (30)  in  nondimensional  form  are 
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WIBWWHIUJ* 


p  * 


Y(Y-i)M2 


T 


P 

P 


It  can  be  seen  that  the  general  viscous  hydromagnetic  equations 
contain  six  fundamental  nondime  ns  ional  parameters.  These  are 


M 

S 


m 


R 


Mach  number 
Interaction  parameter 


u0/<vp0/p0)1/2 

(B^/2®0)/(p0U^) 


Magnetic  Reynolds  number  (UqL/tJq) 
Viscous  Reynolds  number  PqUqL/uq 


R 


Viscous  Prandtl  number  (C  UqAq) 


We  note  that  the  Alfv£n  speed  is  defined  as 

1/2 

CA  ■  (^V) 

and  the  Alfven  Mach  number  is  M^  =  U/C^.  Hence  the  interaction  param¬ 
eter  S  is  also  twice  the  reciprocal  of  the  square  of  the  Alfven  Mach  number. 

For  flows  in  the  absence  of  viscous  and  diffusion  effects 
Rm  “*■00,  Rg  —  oo  we  obtain  the  inviscid  hydromagnetic  equations: 


i£  +  V  *  M  =0  <47) 

at 


tiL  +  v  •  [  5”  *  s  T*  1  =  0  (48) 

at  L  J 


12 


8e 

at 


+ 


(49) 


—  ■  VXfUXB’rO  (50) 

at 


The  jump  equations  across  hydromagnetic  shocks  immediately  follow 
from  Eqs.  (47) -(50). 


(51) 


(52) 


(53) 


(54) 


where  {  }  denotes  the  difference  in  the  quantity  across  the  shock  surface 
and  IT  is  the  normal  to  the  shock  surface. 

Since  the  electrical  conductivity  achievable  in  nonideal  plasma 
is  large  but  finite,  the  magnetic  Reynolds  numbers  are  not  infinite  but 
perhaps  vary  in  the  range  1  S  Rm  S.  20. 

In  this  range,  the  appropriate  system  is  the  inviscid,  finite 
conductivity  hydromagnetic  system: 

+  V  ■  M  =0  (55) 

at 

£2-  +  V  •  f  =0  (56) 

at 
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or  in  the  alternative  form,  the  momentum  and  energy  equations  are 


2#  ♦  V  •  5  .  s  7  •  r, 

dt 


9e 

at 


— X  ^  «  i  ^  <v 

+  v  •  h  =  s  v  •  (U  •  r*)  -  SRmv  •  ( n  v  •  n 


(56a) 

(56b) 


2.  6  Applied  and  Induced  Fields 

Let  us  separate  the  magnetic  field  B  into  an  applied  portion 
B  ^  sustained  by  currents  external  to  the  plasma  and  plasma  induced 
portion  B  ^  which  results  from  currents  flowing  within  the  plasma: 

B  =  B(0)  +  B  (i)  (57) 


The  induction  equation,  Eq.  (58),  then  becomes 

a TT(i)  —  17  n(i)  i  f — 

2%--Vx(UXB  =  -R  *Jvx 

at  m  l 


^  ^  *w  (^)  ^  ^ 

(r?V X B  )  *  Vx(r)  J 


K 


„  T(0)  17(0) 

■R ‘VxfnVxB  )-b 
•  m 
17(0)  ^ 

where  is  denoted  B  ^ 

at 


(58) 


3.0  QUASI-ONE -DIMENSIONAL,  TRANSIENT  MAGNETOGASDYNAMICS 

OF  HYPERVELOCITY  PULSED  FLOWS 

3. 1  Plasma  Flow  Configuration 

We  now  consider  the  behavior  of  shock  generated  magneto¬ 
hydrodynamic  interaction  and  low  to  high  magnetic  Reynolds  numbers 
according  to  a  quasi-one -dimensional  description.  Such  a  flow  consists 
of  a  hot  plasma  "plasmoid"  formed  between  a  driven  ionizing  shock  wave 
and  its  following  contact  surface.  The  plasmoid  is  created  by  a  sudden 
release  of  energy  in  a  driver  section  which  is  in  contact  with  a  test  gas 
in  which  the  plasmoid  propagates.  Such  a  flow  may  be  driven,  for 
example,  by  the  use  of  focused  chemical  explosives.  [:],  [2] 

The  conducting  plasmoid  enters  a  region  in  which  an  externally 

imposed  magnetic  field  Bq  and  electrodes  coupled  to  an  external  circuit 

exist  (Fig.  3-t).  The  plasma  conducts  current  to  this  external  circuit  and 

is  subject  to  Lorentz  forces  and  Joule  heating  as  it  propagates  through 

the  magnetic  field.  If  the  explosion  drive  is  a  chemical  source,  such  a 

plasmoid  will  be  of  the  order  of  5-20  cm  in  length  in  traversing  a 

magnetic  field  region  of  the  order  of  100  cm  at  velocities  of  the  order  of 
4  -1 

10  ms  .  The  plasmoid  may  exist  at  pressures  up  to  1  k  bar  and 
energies  of  5  eV. 


If  °"q »  Pq *  are  the  characteristic  electrical  conductivity,  mass 
density  and  velocity  within  the  plasmoid,  the  flow  may  be  specified  by  an 
interaction  number  i  and  magnetic  Reynolds  number  rm  (in  addition  to  the 
gasdynamic  Mach  number). 


CONTACT 


Fig.  3* 1 .  Schematic  of  explosion  driven  plasmoid  flow. 


Fig. 3 -2.  Space -time  diagram  of  strong  Interaction  hypervelocity 
plasmoids.  Plasmoid  trajectory  in  the  limit  of  vanishing  interaction 
is  delimited  by  the  shock  front  trajectory  SS' S"  and  contact  surface 
trajectory  C  C1  C".  Plasmoid  shock  front  reflection  S'C'  and  trans¬ 
mission  S'S"  are  Initiated  at  encounter  with  magnetic  field.  Contact 
surface  encounters  reflected  front  at  C*  which  generates  fast  reflected 
shock  II  (C,C,M)»  *low  reflected  shock  II  (C,C,,,,)t  and  convected  wave 
(C«C"). 
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For  an  interaction  region  of  length  L,  the  nondimensional  numbers  are 
defined  as 


,L. 


I  =  f  i  dx 
0 


m 


■/ 


r  dx 
m 


(59) 


When  Rm  >>  1,  the  appropriate  measure  of  the  interaction  is  the 

i 

parameter  S  defined  as 


s' 5  (b^/£0)/<PU2>  (60) 

where  the  spatial  average  (  )  is  over  the  electrically  conducting  portion  of 
the  region  L.  For  a  uniform  plasmoid  of  length  a,  these  numbers  become 

1  =orOB0  a//p0U0*  Rm  =  <r0p0U0a’  S=  ^rn- 

Pulsed  magnetohydrodynamic  flows  have  been  examined  in  the 
case  of  low  magnetic  Reynolds  number  (Rm  <<:  1)  and  weak  interaction 
(I  <1)  [3] ,  [4] .  Because  of  the  low  interaction,  these  studies  revealed  simple 
current  flow  through  the  plasmoid  and  weak  magnetohydrodynamic 
deceleration. 

The  transverse  ionizing  shock-wave  which  forms  the  front  of 
the  plasmoid  has  been  extensively  studied  in  the  limit  of  infinitely  lavge 
magnetic  Reynolds  number  [5] ,  [6] .  In  addition  to  the  exposition  of  the 
general  Rankine -Hugoniot  conditions  for  these  shocks  [5]  it  has  also  been 
demonstrated  that  such  shock  waves  can  be  reflected  as  well  as  transmitted 
upon  encounter  with  an  externally  imposed  magnetic  field.  These  studies 
also  showed  that  the  electric  field  in  front  of  the  shock  must  be  self-consistently 
determined  with  the  dynamical  state  behind  the  shock  and  the  electrical 
boundary  conditions  imposed  upon  the  gas  [6]. 
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In  the  present  study,  we  examine  the  magnetohydrodynamics  of 
the  whole  plasmoid  in  its  encounter  with,  and  transit  through  an  externally 
imposed  magnetic  field.  We  show  that  under  conditions  of  strong  inter¬ 
action,  hypervelocity  plasmoids  can  possess  a  rich  variety  of  magneto  - 
hydrodynamic  phenomena  including  magnetically  reflected  shock  waves, 
embedded  MHD  discontinuities,  and  significant  periods  of  transonic  flow 
within  the  plasmoid.  In  particular,  we  reveal  the  dynamics  of  reflected 
and  transmitted  waves  through  the  plasmoid  in  both  the  low  and  high 
magnetic  Reynolds  number  regime.  We  reveal  the  behavior  of  electro- 
thermally  unstable  plasmoids.  We  show  that,  in  general,  the  plasmoid 
is  not  delimited  by  the  region  between  the  shock  front  and  contact  surface. 
Instead,  the  plasmoid  develops  its  own  internal,  evolving  structure 
governed  by  the  mutual  interaction  of  self-heating  and  self-induced  fields. 

Shock-generated  hypervelocity  flows  of  this  kind  are  subject  to 
a  variety  of  nonideal  phenomena.  These  include  wall  interaction  effects 
(viscous  losses,  gas  leakagei  and  ohmic  voltage  drops  in  boundary  layers), 
thermal  radiation  losses,  and  kinetic/ionization  relaxation  effects  behind  shock 
waves.  In  the  present  study  we  ignore  these  effects  and  examine  those  phenomena 
which  arise  specifically  from  the  magnetohydrodynamic  interaction. 

In  Part  3.2  we  present  the  quasi-one-dimensional  version  of 
the  system  of  equations  discussed  in  Part  2.  In  Fart  3.3  we  examine 
the  dynamics  of  strong  interaction  plasmoids  with  an  applied  magnetic 
field  but  at  low  magnetic  Reynolds  number.  In  Part  3.4  we  similarly 
consider  strong  interaction  plasmoids  but  at  large  magnetic  Reynolds 
number.  In  Part  3,5  we  illustrate  the  behavior  of  "transitional"  plasmoids. 


These  are  flows  in  which  the  plasmoid  enters  the  magnetic  field  at 
relatively  low  values  of  interaction  parameter  and  magnetic  Reynolds 
number.  As  a  result  of  self  Joule  heating,  however,  the  plasmoid  con¬ 
ductance  is  elevated  as  it  progresses  through  the  field  carrying  it  into 
the  strong  interaction,  high  magnetic  Reynolds  number  regime. 

3.  2  The  One -Dimensional  Description 


We  consider  a  quasi-one -dimensional  description  of  the  gas 
moving  over  the  spatial  coordinate  x  in  time  t.  If  the  equations  of 
Section  2  are  averaged  over  the  cross  section  of  the  duct,  we  obtain 
the  quasi-one -dimensional  forms 


aw  a  F. 
at  ~  3x 


(61) 


where  W(x.t)  is  the  vector  of  mass,  momentum,  and  total  energy  densities 


W  (x,  t) 


P 

m 

e 


(62) 


In  the  above,  the  following  definitions  apply 

m  =  pU 
e  =  p(  t  +  U2/2 

where  all  quantities  are  to  be  interpreted  as  averages  over  the  duct 
cross  section.  The  magnetic  field  is  that  induced  by  the  plasma 
currents  where  Bq  is  the  applied  field: 

B  =  Bq  +  Bi 
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The  convected  fluxes  of  mass,  momentum,  and  energy  are  contained  in  the 
vector  F  while  H  contains  the  Lorentz  force  and  power  associated  with 
the  magnetic  field  and  the  Joule  dissipation.  These  are  expressed  in 
terms  of  current  density  J  and  magnetic  field  B  : 


m 

0 

F  = 

m2/p  +  P  +  B2/  2{xq 

<TxVx 

m/p(e  +  p) 

J2/cr  -  J  •  (UX  B  ) 

In  this  illustration  study  we  assume  that  the  kinetic  effects  are  confined 
to  the  relaxation  layer  at  the  shock,  and  further,  that  the  relaxation 
layer  is  thin  compared  to  the  overall  thickness  of  the  plasmoid. 

For  the  simple  geometry  (x,y,z)  of  Fig. 3-1,  the  magnetic  field 
B  is  given  by  B(0,0,B),  the  electric  field  by  E  =  E(0,E,0),  and  the 
current  by  J  =  J(0,J,0).  The  description  for  the  near  fields  J,B^  of  the 
plasmoid  is  then  given  by  Ohm's  law  and  the  Maxwell  equation  Ln  the 
MHD  approximation: 


J  =  or  (e  -  UB.  -  UBq)  (64) 

3B. 

*3“  =  -n( E  -  UB.  -  UB0)  (65) 

where  ri  =  (p^o-)  *  is  magnetic  diffusivity. 

External  interaction  conditions  with  an  external  circuit  including 


inductive  coupling  with  the  applied  magnetic  field  coil  are  required  to 


complete  the  description  of  actual  flow  situations.  Rather  than  include 
such  circuit  detail  in  these  illustrations  we  assume  that  the  external 
circuit  is  configured  so  that  an  electric  field  E  =  E(0,E,0)  is  maintained 
within  the  interelectrode  region  whose  magnitude  is  uniform  in  space  and 
given  by 

E  =  /c(  UB) 

where  K  is  a  "load"  parameter  (OS  K  <  1).  For  a  passive  external 
circuit,  the  value  K  =  0  corresponds  to  a  shorted  external  circuit;  for 
K  =  1  the  external  circuit  is  open  circuited. 

The  fluid  variables  p,p,  T,<r,U  are  nondimensionalized  by  the 
values  pQ,  Pq,  Tq,  o-q,  Uq  characterizing  the  interior  of  the  initial  plasmoid 
before  encounter  with  the  magnetic  field.  Nondimensional  space  and  time 
x,T  are  defined  in  terms  of  x,  nondimensionalized  by  the  plasmoid  length 
a,  and  t  by  a/Ug.  The  nondimensional  parameters  governing  the  inter¬ 
action  are  the  Mach  number  M,  the  gas  heat  ratio  y,  either  of  the  inter- 

I 

action  parameters  I  or  S,  and  the  magnetic  Reynolds  number  Rm. 

Boundary  and  Initial  Conditions 

In  the  limit  of  Rm  — •  oo,  the  system  consisting  of  Eqs.59,  and 
(62)  is  fully  hyperbolic.  For  finite  R  ,  the  system  is  mixed  hyperbolic/ 
parabolic  with  embedded  regions  where  resistive  effects  occur.  The 
boundary  and  initial  conditions  which  specify  the  interaction  problem  for 
an  explosion  generated  plasmoid  encounter  with  a  magnetic  field  are  as 
follows.  As  an  initial  condition  we  take  an  idealized  explosion  driven 
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flow  in  which  the  plasmoid  of  given  breadth  a  occupies  the  hot  zone 
between  contact  surface  and  shock  front  [  9  ]  (Fig.  3-3).  At  time  t  =  0 
the  shock  front  is  located  at  the  edge  of  the  magnetic  field.  Over  the 
time  scale  for  the  dynamics  of  interest  the  shock  front  of  the  plasmoid 
will  run  continuously  into  the  quiescent  driver  gas  while  the  backward 
running  rarefaction  continuously  runs  into  the  explosive  source.  Hence 
the  boundary  conditions  for  the  fluid  equations  are  those  of  specified 
explosion  and  quiescent  states  at  the  boundaries  x  =  +L^,  x  =  - L ^ 
respectively. 

The  boundary  condition  for  the  induced  magnetic  field  from 
Eq.  (61)  is  that  of  symmetry  across  the  overall  plasmoid  so  that  atx  =  L^, 
x  =  -L^  which  lie  outside  the  region  of  any  current  flow 

B.(-L2,t)  =  -B.O^t) 

The  applied  magnetic  field  Bq  is  uniform  in  both  space  and  time. 

Numerical  Procedures 

The  solutions  to  the  initial-value  problems  formulated  above  and 

to  be  discussed  in  Section  3.3,  3.4,  and  3.5  are  computationally  generated 

with  second  order  accurate  explicit  finite  difference  operators.  The 

hyperbolic  system  is  treated  with  the  MacCormack  version  of  the  Lax- 

Wendroff-Richtmyer  operator  [  10  ] .  For  the  space-time  grid  utilized, 

comparisons  were  made  with  the  analytically  available  solutions  for  the  zero 

and  infinite  R _ _  zero  interaction  limits.  At  the  extreme  pressure  ratios  of 

m 

5 

10  between  driver  and  driven  gas  for  these  explosion  generated  plasmoids. 


the  maximum  variations  between  the  computationally  generated  and  analytical 
solutions  within  the  plasmoid  (expressed  as  a  fraction  of  the  analytical  solution) 
are  0.025  in  velocity,  0.04  in  pressure,  and  0.08  in  temperature. 


STD  9-2976A 


Fig.  3-3.  Initial  condition  for  plasmoid  interaction  for  the  pressure  (a), 
temperature  (b),  and  velocity  (c).  Shock  front  of  plasmoid  of 
given  breadth  a  located  at  magnetic  field  edge  at  time  t  =  0 

STD  RESEARCH  CORPORATION 
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Interaction  at  Low  Magnetic  Reynolds  Number 


We  now  proceed  to  the  first  of  several  illustrations  of  the  fore¬ 
going  description.  We  consider  first  the  strong  interaction  of  a  plasmoid 
with  the  magnetic  field  but  at  low  magnetic  Reynolds  number.  In  Fig.  3-2 
the  kinematics  of  this  situation  are  shown.  When  the  incident  plasmoid 
encounters  the  magnetic  field,  the  leading  shock  front  may  be  both  trans¬ 
mitted  and  reflected.  In  the  case  of  reflected  fronts,  the  rear  (contact 
surface)  of  the  plasmoid  subsequently  interacts  with  the  reflected  shock 
front.  This  colliding  disturbance  then  radiates  a  fast  and  slow  reflected 
shock  (d.-noted  shock  II)  back  through  the  plasmoid  (which  consists  of 
subsonic  flow  behind  the  reflected  shock  front)  where  it  then  collides  with 
the  now  strongly  decelerated  shock  front. 

For  the  illustration  shown  here,  we  select  a  plasmoid  with  Mach 
number  M  =  1.64,  interaction  parameters  1  =  20,  S  =  200,  and  Reynolds 
number  Rm=0.1,  just  before  encountering  the  magnetic  field.  The  full 
conditions  for  the  flow  are  given  in  Table  I.  We  impose  the  condition 
that  the  electrical  conductivity  is  spatially  uniform  within  the  high  tempera¬ 
ture  plasmoid  (we  consider  electrical  conductivity  functions  which  are  con¬ 
sistently  coupled  to  the  gas  thermodynamic  state  in  Part  3.4).  This  uniform 
conductivity  distribution  is  achieved  in  the  computations  with  the  model 
conductivity  function 


T/T0  <  1/3 


T/Tq  >  1/3 
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TABLE  I 


Conditions  for  Interaction  at  Low  Magnetic  Reynolds  Number 
T  =  Quiescent  driven  gas  temperature 

CO 

2  Quiescent  driven  gas  pressure 


M  =  1.64 

V  = 

1.5 

Po/Poo  =  232 

I  = 

R  = 
m 

20 
0.  1 

Tft/T  =  45 

O'  00 

K  = 

0.5 

STD  9-3050 


TABLE  II 


Conditions  for  Interaction  at  High  Magnetic  Reynolds  Number 


M  =  1.64 
P0/Poo  =  232 

VTco  =  45 


y  =  1.5 

I  =  20 

=  5 
m 

K  =  0.5 


STD  9-3051 


TABLE  III 


Conditions  for  Transitional  Plasmoid 


M  *  1.64 

P0/P«  *  232 

VTcc  *  45 

V  =  1.5 


<rQ  =  1080 
n  =  3.11 
Z"2£f/fnA  =  0.16 
K  =  0.5 


STD  9-3052 


which  effectively  switches  on  a  constant  conductivity  <Tq  within  the  plasmoid  and 
switches  the  conductivity  off  outside  the  zone  in  which  the  plasmoid  exists.  The 
dynamics  of  the  interaction  are  exhibited  in  Figs.  3-4  and  3-5.  When  the 
plasmoid  shock  front  encounters  the  magnetic  field,  it  is  both  reflected  and 
transmitted.  The  reflected  shock  I  collides  with  the  contact  surface  and  initiates 
a  fast  reflected  shock  II  and  a  slow  reflected  shock.  The  fast  reflected  shock  II 
reestablishes  high  velocity  flow  through  the  plasmoid  and  reencounters  the 
transmitted  shock  front  which  has  been  decelerated.  During  this  period  of 
strong  wave  dynamics  the  current  distribution  within  the  plasmoid  is  strongly 
affected  (Fig.  3-5).  The  current  is  diminished  to  very  small  values  during 
the  period  of  plasmoid  deceleration  behind  the  reflected  shock  I,  and  the 
returns  back  to  enhanced  levels  after  passage  of  the  fast  reflected  shock  II. 

The  low  magnetic  Reynolds  number  of  the  plasmoid  allows  the  current  to 
diffuse  nearly  uniformly  throughout  its  breadth. 

3. 4  Interaction  at  High  Magnetic  Reynolds  Number 

We  next  consider  the  behavior  of  a  uniform  conductivity  plasmoid 
in  the  high  magnetic  Reynolds  number  regime.  For  this  case  the  incident 
plasmoid  has  an  interaction  parameter  I  =  20  as  in  the  previous  illustration 
but  a  magnetic  Reynolds  number  Rm  =  5.  Correspondingly,  the  interaction 
number  S'  has  the  reduced  value  S'  =4.  The  full  conditions  for  this  flow 
are  given  in  Table  II.  The  encounter  of  this  plasmoid  with  the  magnetic 
field  is  similar  to  that  of  Part  III.  The  features  peculiar  to  the  higher 
Reynolds  number  are  best  perceived  in  the  current  distribution  of  Fig.  3-7 
which  is  more  nonuniform  compared  to  that  of  Fig.  3-8.  When  the 
current  levels  rise  behind  the  reflected  shock  II,  they  do  so  by  directly 
following  the  shock  until  it  merges  with  the  transmitted  front.  The  current 
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(c)  REFLECTED  SHOCK  II  TRAJECTORY 
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b)  PLASMOIO  INITIALLY  HEATED  AT  LOU  CURRENT 
LEVELS  AND  CONVECTED  AT  LOW  VELOCITY 
BEHIND  REFLECTED  SHOCK  FRONT  I 
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(a)  PEAK  IN  CURRENT  INITIATED  AT 

ENCOUNTER  OF  PLASNOIO  SHOCK  FRONT 
WITH  MAGNETIC  FIELD 
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absent)  behind  shock  front 


(b)  CONTACT  SURFACE  ENCOUNTER 

WITH  MASNETIC  FIELD  INITIATES 
A  SLOW  WAVE,  FAST  WAVE,  A NO 
CONVECTE0  WAVE  v 


ce**  i-s 

?  O  >m  D  J 

fi  b  O  J3  I 

■M  c  ^  ! 

*  01  01  C  -w  s 

— *  ,  01  3  V*  I 

»  V  5  c  §  J 

wfo»r; 

t  “g §  3 ; 

2  ^  aoiJ  T3  1 

<n  c  Eh  Ka  < 

^  W  L  D  C  ( 
0>  *  JCS  3 

rTJ  *! 

4-»  fu  fli  1 

*^o«" 

C  Q  £j  tw  > 

bOTJ  £  A  V  J 

«  »  M0  -S'- 

£  o  «  .-2 

2  i  «T  5 

2  <u  *?  ^ 

■W  b<  5  O  *r!  I 

«  ouS  a  i 

q-J  C  ^4  fli  O  ^ 

2  *  IJ  1 

°  «  a,2  m  ‘ 

H  TJ  c  01  &  i 

9  O  O  “ 
(0  oi  tj  ai  "j  -j 
n  15  o  -o  53  .1 
^  «  « 


! 


j 


31 


maximum  then  remains  at  the  shock  front  of  the  plasmoid.  As  a  result  o 
decelerating  Lorentz  forces  concentrated  immediately  behind  the  shock 
front,  the  shock  front  is  slowed  and  the  overall  breadth  of  the  plasmoid 
is  decreased  as  it  progresses  through  the  magnetic  field.  This  is  in 
contrast  to  the  plasmoid  dynamics  of  Parts  3.3  and  3.5. 

3. 5  Transitional  Plasmoids 

We  now  turn  to  consideration  of  plasmoid  behavior  with  a 

coupled  electrical  conductivity  model.  In  contrast  to  the  previous 

illustrations  in  which  the  conductivity  is  spatially  uniform  within  the  high 

temperature  plasmoid  and  vanishes  outside,  we  consider  a  conductivity 

which  is  appropriately  coupled  to  the  thermodynamic  state  of  the  gas. 

As  a  result,  local  regions  within  the  plasmoid  can  be  rendered  more 

conductive  by  the  self  Joule  heating  of  the  plasmoid.  With  the  electrical 

conductivity  strongly  coupled  to  the  Joule  dissipation,  a  plasmoid  in  the 

low  I,  low  R  range  can  evolve  into  the  large  I,  large  R  range  as  it 
m  m 

progresses  through  the  magnetic  field  and  experiences  further  Joule 
heating.  We  term  such  flows  "transitional"  plasmoids. 

We  consider  a  conductivity  function  of  the  form 


-1 

(T 


In  the  above  <r  is  an  electrical  conductivity  of  a  neutral  species  back- 
en 

ground  and  <r  .  is  the  Coulomb  conductivity.  We  use  as  a  summary 
representation  of  these  two  contributions  the  forms 
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<r  .  =  152  2  „  T^/fnA 
ei  eff 

where  is  the  average  effective  ionic  charge  and  o-Q.Tg.inA,  n>0  are 

parameters  for  a  given  gas. 

The  conductivity  function  Eqs .  (63)-(64)  is  dominated  by  the 
partially  ionized  conductivity  <ren  at  low  temperature  and  goes  over  to  the 
Coulomb  (fully-ionized)  conductivity  at  high  temperature.  This  conductivity 
function  has  the  property  3<r/3T  e  0  over  the  entire  range  of  temperature. 
Since  there  are  no  thermal  energy  loss  mechanisms  (which  would  be 
principally  radiative)  included  in  the  model,  the  plasmoid  is  unconditionally 
electro  thermally  unstable  [ll],[lz].  This  convective  instability  is  simply  a 
growth  of  temperature  nonuniformities  within  the  plasmoid  due  to  intensifying 
Joule  heating  resulting  from  growing  electrical  conductivity. 

The  interaction  of  a  representative  transitional  plasmoid  is  shown  in 
Figs.  3-8  to  3-10.  This  plasmoid  has  interaction  parameters  1=1,  S  =  1  and 
magnetic  Reynolds  number  Rm  =  1  just  before  it  enters  the  magnetic  field. 
The  complete  conditions  .for  this  flow  are  given  in  Table  III.  It  should  be 
noted  that  the  interaction  at  entry  into  the  magnetic  field  is  considerably 
smaller  than  the  interaction  described  in  Farts  3.3  and  3.4.  The  plasmoid 
progresses  into  the  field  where  it  begins  to  self-heat  and  decelerate.  The 
modest  interaction  at  plasmoid  entry  to  the  magnetic  field  does  not  create 
distinct  reflected  waves,  but  rather  a  strong  and  continuous  deceleration. 
Magnetic  Reynolds  number  and  interaction  parameter  grow  significantly  as 
the  plasma  is  heated.  At  time  ~  =  2. 6,  the  magnetic  Reynolds  number  is 
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(c)  ENCOUNTER  OF  ELECTROTHERMAL 
INSTABILITY  WITH  FAST  WAVE 


-1C.  Temperature  field  corresponding  to  Figs.  3-8,  3-9  showing  developing  electrothermal  instability 


k 


in  excess  of  10  and  the  current  has  progressively  become  sheet-like.  It 

should  be  noted  that  the  current  maxima  no  longer  follow  immediately 
behind  the  transmitted  shock  front.  Rather,  the  current  concentrates  in 
the  electrothermally  heated  zone  which  is  then  convected  at  the  local  fluid 
speed  rather  than  radiated  at  the  shock  front  speed.  This  decelerating 
electrothermal  instability  is  then  swept  up  by  collision  with  the  waves 
initiated  by  the  arrival  of  the  contact  surface  at  the  magnetic  field  inlet. 

A  feature  of  note  is  the  development  of  reversed  current  flow  in  the  up¬ 
stream  portion  of  the  current  sheet  due  to  the  sharply  diminished  magnetic 
field  behind  the  current  sheet  at  large  R  . 

3.  6  Summary  Remarks 

In  this  study  we  have  illustrated  significant  purely  magnetogas  - 
dynamic  phenomena  which  occur  when  a  hypervelocity  pulse  of  plasma 
("plasmoid")  encounters  an  applied  magnetic  field  under  strong  interaction 
conditions.  With  uniform  electrical  conductivity  within  the  plasmoid  (and 
vanishing  electrical  conductivity  outside),  reflection  and  transmission  of 
the  plasmoid  shock  front  are  possible  coupled  with  strongly  nonuniform 
current  evolution  in  time.  Such  plasmoid s  with  large  magnetic  Reynolds 
numbers  have  current  distributions  (and  decelerating  Lorentz  forces) 
concentrated  immediately  behind  the  shock  front.  As  a  result  of  shock 
front  deceleration,  these  plasmoids  diminish  in  breadth  as  they  proceed 
through  the  magnetic  field.  With  electrical  conductivity  within  the  plasmoid 
coupled  to  its  thermodynamic  state,  the  plasmoid  is  electrothermally  un¬ 
stable  and  creates  its  own,  evolving  region  of  enhanced  electrical  conductance 
which  carries  most  of  the  current  and  is  convected  at  the  fluid  speed  within 


r 

I 


i 


I 


the  plasmoid.  The  shock  front  becomes  increasingly  free  of  current  and 
runs  progressively  farther  ahead  of  the  unstable  current  structure  embedded 
in  the  plasmoid  interior. 

Nonideal  phenomena  such  as  viscous  wall  layers,  kinetic  - 
relaxation  effects  behind  the  shock  front,  and  thermal  radiation  losses 
can  play  important  roles  in  these  strong  interaction  flows.  The  basic 
structure  of  the  magnetohydrodynamic  interaction  itself,  however,  is  a 
prerequisite  to  the  description  and  understanding  of  these  additional 
modifying  effects. 


| 
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4.  CURRENT  AND  MAGNETIC  FIELDS  IN  TWO-DIMENSIONAL 


4.  i  Channel  and  Applied  Magnetic  Field  Configuration 

We  now  consider  some  illustrative  flows  in  which  the  fluid 
distribution  of  velocity  and  electrical  conductivity  are  specified  as 
functions  of  space  and  the  induced  magnetic  fields  and  the  plasma 
currents  are  to  be  determined.  Three  general  classes  of  flow  will  be 
examined . 

(a)  Uniform  velocity  and  conductivity  distributions 

(b)  Nonuniform  velocity  and  conductivity  distributions  resulting 
from  an  hypothesized  oblique  shock  system  within  the  MHD 
generator  duct 

(c)  Nonuniform  velocity  and  conductivity  distributions  resulting 
from  supersonic  boundary  layers  on  the  walls  of  the  MHD 
generator  duct. 

All  flows  exist  within  the  duct  geometry  and  applied  magnetic 
field  distribution  shown  in  Figs.  4-1  through  4-10  in  which  the  generator 
electrode  length  is  equal  to  the  duct  height.  The  magnetic  field  is  given 
by 


where  is  the  maximum  applied  field.  The  values  selected  are 

m 

L  =  4h  and  a  =  2/h. 


Since  these  calculations  decouple  the  electricity  from  the  fluid 

behavior  (weak  or  vanishing  interaction)  only  the  magnetic  Reynolds 

number.  R  .  is  a  relevant  electrical  parameter.  All  cases  include  con- 
m 

ditions  for  a  magnetic  Reynolds  number  Rm  >  1  which  is  the  range  of  interest. 
4.  2  Two-Dimensional  Electrical  Description 


Let  us  now  consider  the  forms  of  Eq.  (58)  appropriate  to 
conduction  in  the  plane  perpendicular  to  an  applied  magnetic  field 
Let  the  z  axis  be  aligned  with  the  applied  field  B  ^  and  let  x,  y  be  the 
coordinates  defining  the  plane  of  conduction.  The  current  vector  then 

becomes  J  =  J(J  . J  ),  the  magnetic  vector  becomes  B  =  B(0,0,B),  and 
x  y 

the  velocity  U  =  U  (U^,  U^.,0),  with  Eq.(iO)  caking  the  form 

*»  ag(i)  ao(i) 

Jx  =  J  =  -  %-  (66) 

X  9x  y  9y 

The  induction  equation  (58)  is  then 


9B 


«(i> 


9t 


+  4(UxB(i|)  +  ~(UvB(i))  =  R 
9x  X  9y  y  m 


*  L11  ( 


2b1  92B^\ 

9x2  +  By2  ' 


*  4^’  ^^-'1-4  <UB«»>  -  j-  (U  B 
9x  9x  9y  9y  -I  9x  Y  Y 


)rj  dB(i) 

By  9y 


5(0) 


) 


Bq  (67) 


The  boundary  conditions  appropriate  to  Eq.  (67)  are  that  on 
insulating  boundaries. 


constant 


(68) 


while  on  conductors 


where  n  is  the  normal  coordinate  to  the  conductor  surface. 


(69) 
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At  steady-state  the  electric  field  is  derivable  from  a  potential  #: 


E  =  -V* 


(70) 


5(i) 


Given  the  induced  fields,  B  '  ,  the  potential  distribution  $(x)  is 


®(x)  =  *(x0)  +  £  ['Rm(^X  ®(i))  +  UX  B  j  •  d 


(71) 


-x  r 


where  d  s  is  the  line  integral  along  any  contour  from  Xq  to  x. 


4. 3  Uniform  Velocity  and  Electrical  Conductivity 

4.3.1  Load  with  Point  Electrode 


If  the  flow  discussed  in  4.3.1  is  suggested  with  point  electrodes 
and  a  load  current  per  unit  depth  I  is  passed  through  the  circuit,  the 
field  distribution  shown  in  Fig.  4-1  results.  These  results  for  B^/B^ 
are  for  a  current  I  s  B^'  =  0.25  and  a  Reynolds  number  RM  =  10. 

The  convection  of  the  magnetic  field  down  stream  by  the  fluid  is  con¬ 
siderable  and  the  eddy  current  cells  at  the  magnetic  field  edge  are  as 
significant  as  the  generator  current.  The  result  of  Gill  [13]  is  con¬ 
sistent  with  this  computation. 


4.3.2  Load  with  Finite  Electrode 


If  the  basic  situation  above  is  supplied  with  finite  electrodes, 
the  upstream  eddy  current  cells  couple  its  current  into  the  generator 
circuit  as  shown  in  Fig.  4-2.  If  the  generator  current  is  increased  to 
1=1,  the  result  is  shown  in  Fig.  4-3. 
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4.4  Nonuniform  Velocity  and  Conductivity  Distributions  Resulting 

from  a  Shock  System 

We  now  consider  nonuniformities  in  electrical  conductivity  and 
velocity  resulting  from  an  hypothesized  oblique  shock  system.  We 
consider  discontinuous  distribution  across  the  shock  with  (  ) ^  denoting 
the  upstream  side  of  the  shock  and  (  ) ^  denoting  the  downstream  side. 

The  magnetic  Reynolds  number  based  upon  upstream  conditions  is  R. .  -  10 
for  all  cases.  The  current  is  I  =  1  for  all  cases. 

4.4.1  Shock  in  front  of  Electrodes 


When  the  shock  is  in  front  of  the  electrodes,  the  induced  field 
distribution,  shown  in  Figs.  4-4  through  4-6  result.  In  Fig.  4-4  there 
is  a  conductivity  group  of  3  and  no  velocity  jump.  In  Fig.  4-5  the 
conductivity  jump  is  10  with  no  velocity  jump.  In  Fig.  4-6,  the  con¬ 
ductivity  jump  is  10  and  there  is  a  velocity  jump  of  i/2. 

4.4.2  Shock  in  Channel  Center 

In  Fig.  4-7,  the  conductivity  jump  of  10  is  shown  with  the 
shock  system  located  in  the  channel  center.  With  the  shock  further 
downstream,  the  front  eddy  cell  is  less  free  to  couple  into  the  electrodes 
and  generator  circuit. 

4.4.3  Shock  at  Upstream  Edge  of  Magnetic  Field 

When  the  shock  is  moved  upstream  to  the  upstream  edge  of 
the  magnetic  field,  the  electrodes  strongly  couple  the  upstream  eddy  cell. 

4. 5  Nonuniform  Velocity  and  Conductivity  Distributions  Resulting 

from  Supersonic  Boundary  Lavers  on  the  Walls  of  the  MHD  Generator 
Duct 

4.5.1  Flow  Structure  for  Electrical  Interaction 


The  flow  field  has  been  calculated  as  that  resulting  from  a 
reservoir  which  feeds  the  duct  with  Argon  at  an  inlet  velocity  of  10  km/sec, 
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an  internal  energy  of  37  Mj/kg,  a  pressure  of  10  k  bar,  and  a  nominal 
electrical  conductivity  of  25,000  mho/m  at  the  duct  inlet.  The  resulting 
turbulent  boundary  layer  interaction  is  considerable  as  the  flow  proceeds 
down  the  duct. 


The  generator  electrodes  are  located  at  0.6  m  downstream 
from  the  duct  entry  from  the  driver.  At  the  station  10  cm  upstream 
from  the  generator  inlet  (x  =  0.50  m)  the  conditions  of  the  flow  are 


Velocity  at  duct  centerline 
Temperature  at  duct  centerline 
Mach  number 
Boundary  layer  thickness 


9612  m/sec 
36,800  K 
3.14 
5. 5  mm 


At  the  center  of  the  generator  section,  (x  =  0.60  m)  these  values  are 


Velocity  at  duct  centerline 

9520  m/sec 

Temperature  at  duct  centerline 

36,930  K 

Mach  number 

3.11 

Boundary  layer  thickness 

6.  5  mm 

.t  of  the  duct,  (x  =  0.70  m)  these 

values  are 

Velocity  at  duct  centerline 

9430  m/sec 

Temperature  at  duct  centerline 

36,970  K 

Mach  number 

3.07 

Boundary  layer  thickness 

7 . 5  mm 

This  distribution  of  nonuniform  velocity  and  temperature  in  x  and  y  (with 
corresponding  density  and  conductivity  nonuniformity)  resulting  from  the 
boundary  layers  is  used  as  the  basis  for  the  electrical  calculations 
utilizing  Eq.  (67).  The  nondimensionalized  values  are  based  upon  the 
velocity  and  conductivity  at  the  duct  inlet  (x  =0).  The  Reynolds  number 
is  based  upon  the  channel  height,  h  and  has  the  value  7.1. 
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4.5.2  Electrical  Conduction  at  Vanishing  Magnetic  Reynolds  Number 

We  first  exhibit  the  nonuniform  electrical  conduction  in  the  channel 
at  low  magnetic  Reynolds  number,  Rj^  =  0.01.  The  total  generator 
current  per  unit  depth,  I,  has  the  value  0.01. 

~~(i) 

The  induced  magnetic  field,  B '  nondimensionalized  on  the 
maximum  value  of  the  applied  field  is  shown  in  Fig.  4-9  corresponding 
to  the  current  I.  It  can  be  seen  there  is  negligible  convection  of  the 
current  distribution;  the  current  spreads  out  to  fill  the  central  portion 
of  the  channel.  The  principal  effects  of  the  electrical  nonuniformities 
are  the  voltage  drops  through  the  cool  boundary  layers  and  the  fringing 
of  the  current  at  the  electrode  edges. 

4.5.3  Electrical  Conduction  at  High  Magnetic  Reynolds  Number 

When  the  magnetic  Reynolds  number  is  R  =7.13  the  resulting 
»  m 

induced  magnetic  field,  B  ,  is  shown  in  Fig.  4-10.  For  the  case  of 
I  =  1  it  can  be  seen  that  the  generator  current  is  driven  downstream  in 
the  usual  fashion,  the  bulk  of  the  power  actually  being  produced  some¬ 
what  downstream  of  the  electrodes.  Because  of  the  cool,  poorly  con¬ 
ducting  boundary  layers,  the  eddy-current  cells  at  the  magnet  edges  do 
not  couple  into  the  generator  circuit  as  they  do  in  the  cases  of  Figs.  4-3 
through  4-8  which  do  not  include  boundary  layer  effects. 
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